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By means of atomistic tight-binding calculations, we investigate the transport properties of ver¬ 
tical devices made of two incommensurately misoriented graphene layers. With a chosen transport 
direction (Ox-axis), we define two classes of rotated graphene lattice distinguished by the different 
properties of their lattice symmetry and, hence, Brillouin zone, i.e., the two Dirac cones are located 
either at the same /c^-point {Ky — Ky = 0) or at different /c^y-points (Ky = —Ky = 27t f^Ly, where 
Ly is the periodic length along the Oy axis). As a consequence, a misalignment of Dirac cones of two 
layers occurs and a significant energy-gap (~ a few hundreds of meV) of transmission is achieved in 
devices made of two layers of different lattice classes. We also shown that strain engineering can be 
used to strongly enlarge the gap in this type of device. 

PACS numbers: xx.xx.xx, yy.yy.yy, zz.zz.zz 


Graphene is nowadays one of the most attractive ma¬ 
terials in several research fields because of its unusual, 
and in many respects, excellent physical properties, as 
a consequence of its two-dimensional honeycomb lattice. 
In particular, it is expected to be a good channel mate¬ 
rial in high-frequency electronic devices, flexible devices, 
spintronic devices, and to provide outstanding proper¬ 
ties for photonics and optoelectronics, sensors, energy 
storage and conversion and so on [T]. However, with a 
view to digital applications in electronics, the lack of 
bandgap in graphene is known to make it difficult to 
turn off the current, which leads to low Ion/Iqff ra¬ 
tio and poor current saturation 0 0, which are im¬ 
portant limitations for innovative applications of tran¬ 
sistors. So far, many efforts of bandgap engineering in 
graphene have been made to solve these issues. For in¬ 
stance, techniques such as cutting 2D graphene sheet into 
narrow nanoribbons [4], depositing graphene on hexago¬ 
nal boron nitride substrate [5], nitrogen-doped graphene 
[6], applying an electric field perpendicularly to Bernal- 
stacking bilayer graphene [7], graphene nanomeshes [8], 
using hybrid graphene/hexagonal boron-nitride [9] or 
vertical graphene channels m have been explored. Al¬ 
though they are certainly promising options for opening 
a bandgap in graphene, some of them still have their 
own issues while the others still need experimental ver¬ 
ification and realization. Hence, bandgap engineering is 
still a timely topic for the development of graphene in 
nanoelectronics. 

In this regard, we have recently proposed the use of 
strain engineering to open a finite energy gap in graphene 
hetero-channels USHIIl- It is basically due to the fact 
that although it can not change the gapless character, 
a small strain can lead to a significant displacement of 
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Dirac cones of the graphene’s bandstructure in the k- 
space. Due to this effect, the misalignment of Dirac cones 
in different graphene sections of a heterochannel occurs 
and may result in a finite energy-gap of the transmission. 
We have hence demonstrated that strain engineering is 
indeed a promising option to enlarge the applications of 
graphene, e.g., as transistors or strain and thermal sen¬ 
sors. 

More interestingly, devices made of vertical stack of 
misoriented (twisted) graphene layers have been inves¬ 
tigated in m- As suggested, this design strategy 
may have the advantages to be based on the use of 
graphene materials only and of making it possible to ap¬ 
ply a uniform strain over the full device channel, com- 



FIG. 1: Schematic of devices made of vertical stack of mis¬ 
oriented graphene layers. The left layer has a zigzag (a) or 
amrchair orientation (b) along the transport direction. 
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FIG. 2: (E — /cy)-maps of transmission probability around the neutrality (Dirac) point in different devices. The devices made of 
(AM 7 , 4 @AM 4 , 7 ), (AM 6 , 6 @AM 7 , 4 ), (ZZ 4 , 4 @ZZ 5 , 3 ), (AMi 3 , 7 @AMi 3 , 8 ) laycrs are considered in the panels (a,b,c,d), respectively 
(see the definition of layers in the text). 


pared to the vertical devices m and strain hetero¬ 
channels na ig. However, we have previously lim¬ 
ited our investigation to devices made of commensu- 
rately twisted layers, an ideal case of twisted multilayer 
system. This naturally gives rise to a question about 
the transport properties of devices made of incommen¬ 
surate graphene layer stacks. Hence, in the current 
work, we aim at answering this question and explore a 
new possibility of energy-gap engineering. Experimen¬ 
tally, controlling the misorientation of graphene layers 
has been explored and used to achieve resonant tun¬ 
nelling in vertical devices based on graphene/hexagonal 
boron-nitride/graphene heterostructures [18]. The com¬ 
mensurate - incommensurate transition in graphene on 
hexagonal boron-nitride has been also achieved in m- 
These works could provide good suggestions for the fab¬ 
rication of the structures investigated here. 

The studied structures are schematized as in Fig. 1. 
They consist in two partially overlaped graphene layers 
and the transport (Ox) direction is perpendicular to the 
overlap section. The top layer is rotated relatively to the 
bottom one by an angle but, differently from the 
study in m, they are incommensurate layers, i.e., they 
do not have any common (truly) periodic lattice vector. 
Here, we consider 2D graphene channels, i.e., the width of 


graphene layers is very large, or, in other word, assumed 
to be infinite. Hence, to carry out the calculations, the 
periodicity of the channel along the transverse (Oy) di¬ 
rection is a necessary assumption [ISHIIl- In this regard, 
we propose an approximation as follows. In all studied 
cases, the two layers have different lengths Ly and Ly 
of unit cells along the Oy axis. When they are coupled 
to form a heterochannel, we can always find two peri¬ 
odic vectors = TUiLiy and = m 2 L^ that satisfy the 
condition Ly Ly. On this basis, our calculations are 
performed, similarly as in [HHn] with an approximately 
periodic length Ly = {Ly + Ly)l2. This approximation 
can be reasonably accepted if the difference between Ly 
and Ly is small enough, as it has been confirmed in simi¬ 
lar systems [TIHI4) . In particular, in studies on graphene 
grain boundaries [HHIS], it was shown that tight bind¬ 
ing calculations made within the above approximation 
are in good agreement with density functional theory, es¬ 
pecially in the region close to the energy-gap. The trans¬ 
port properties have been accurately computed around 
the gap even when the difference between Ly and Ly 
reaches ^ 3.85 %. Similar agreement has been obtained 
in vertical graphene/hBN heterochannels [14], where the 
lattice mismatch is about 1.8 %. Throughout the present 
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FIG. 3: Schematic of two typical unit cells of graphene lattices, shape of their Brillouin zone, and corresponding pictures 
showing their bandstructure prohle along the ky axis. 


work, to ensure the accuracy, our calculations were per¬ 
formed with and Li that have only a small difference, 
i.e., < 2 %. 

Now, we would like to explain the definition of 
graphene lattices considered in this work. In Fig. 1, 
we plot the geometry of two typical devices, i.e., the left 
(bottom) layer has zigzag and armchair orientation in (a) 
and (b), respectively, with respect to the Ox axis. We 
distinguish two lattice types as follows. For both cases, 
the lattice is determined by rotating the original (zigzag 
or armchair) lattice by an angle 0 tl- This angle is com- 

L 

puted as cos 0tl = l|Lylll|L»|| where Ly is the vector (be¬ 
fore the lattice is rotated) defining the size of unit cells 
along the Oy-axis while Ly is the vector of the original 
lattice. These vectors are determined as Ly = giai+g 2 a 2 
and Ly = ai + a 2 . For the type 1 (i.e., ZZqi^q 2 lattice, 
see in Fig. 1(a)) when the original lattice has a zigzag 
orientation along the Ox axis, ai = (—3)ao/2 and 
a 2 = (\/3, 3)ao/2. For the type 2 (i.e., AMgi^g 2 lattice, 
see in Fig. 1(b)) when the original lattice has an armchair 
orientation, ai = (—3, y/3)ao/2 and a 2 = (3,v^)ao/2. 
Here, ao = 0.142 nm is the in-plane C-C bond length 
in graphene. Accordingly, we determine also the size of 
unit cells along the Ox axis as Lx = Piai + P 2^2 that 
satisfies the condition L^, -L^ = 0 , i.e., ^ — and 

^ 2 ;Zgi,g 2 and AM,i,g 2 lattices, respec- 

tively. 

To compute the electronic transport in the consid¬ 
ered systems, we employed atomistic tight-binding cal¬ 
culations and Green’s function techniques as in [m- 
[la MM- In particular, the hopping parameters 
are determined similarly as in [22]-[24], i.e., tij = 
to exip[—P {vij /vq — 1}] where to = —2.7 eV (resp. 0.48 
eF), (3 = 3.37 (resp. 7.42) and ro = 1.42 A = ao (resp. 


3.35 A) for in-plane (resp. interlayer) interactions. 


TABLE I: Detailed description of some devices investigated 
in this work and their transport gap Eg: 


Layer 1 

Layer 2 

4>tl 

Ly (ao) 

7?(%) 

E,{meW) 

ZZ4^4 

ZZ5,3 

8.21° 

12.062 

1.031 

703.2 

ZZs^s 

ZZs,! 

24.18° 

14.899 

1.351 

569.3 

ZZi3,i 

ZZii,4 

11.25° 

23.366 

0.549 

363.0 

ZZ8,8 

ZZi3,2 

22.94° 

24.216 

1.790 

350.2 

ZZi5,l 

ZZi3,4 

9.80° 

26.776 

0.836 

316.7 

ZZi9,i 

ZZi7,4 

7.78° 

33.629 

1.061 

252.2 

AM6,6 

AM7^4 

25.28° 

10.464 

1.369 

810.6 

AMg^g 

AMio,7 

16.99° 

15.491 

1.250 

547.5 

AMi3,7 

AMi3,8 

5.04° 

19.595 

0.781 

432.8 

AMi4^g 

AMi3,11 

12.41° 

21.141 

1.342 

401.2 

AMi4,14 

AMi6,11 

17.78° 

24.402 

1.259 

347.6 

AMi7^9 

AMi7,10 

3.87° 

25.573 

0.458 

331.6 

ZZii,4 

AMi4,13 

18.59° 

23.366 

0.549 

362.9 

ZZi3,2 

AMi6,11 

24.83° 

24.494 

0.500 

346.2 

ZZi5,l 

AMi7,13 

16.20° 

26.776 

0.836 

316.7 

ZZi5,4 

AMi9,14 

26.22° 

29.798 

1.689 

284.6 


First, we would like to analyze some basic trans¬ 
port properties of the considered devices. We display 
in Fig. 2 the {E — /c^)-maps of transmision probabil¬ 
ity in four cases: the devices made of (AM 7 ^ 4 @AM 4 j), 
(AM 6 , 6 @AM 7 ^ 4 ), (ZZ 4 ^ 4 @ZZ 5 ^ 3 ), (AM 13 ^©AMis^g) lay¬ 
ers in the panels (a,b,c,d), respectively. Note that the 
two layers AM 7^4 and AM 4 j are commensurate while the 
other systems are incommensurate. It is clealy shown 
that in the case of commensurate lattices, even if the two 
lattices have different orientations, there is no energy-gap 































4 


of transmission. This can be explained by the fact that 
the commensurate lattices have the same periodic vec¬ 
tors and hence, because of their lattice symmetry, they 
have the same Brillouin zone. Since their Dirac cones are 
located at the same k-point, the transport gap is zero, 
similarly as reported in [nmaiM]. The transport pic¬ 
ture, however, can be dramatically changed in the devices 
made of incommensurate layers, i.e., finite energy-gaps 
can be achieved as shown in Figs. 2(b,c,d). The value of 
energy-gap is quite significant in these cases: Eg SIO.6 
meV in (b), 703.2 meV in (c), and 432.8 meV in (d). This 
phenomenon is explained as follows. For all cases stud¬ 
ied, we can distinguish two classes of rotated graphene 
lattice as shown in Figs. 3(a) and 3(b). In the class 
1 , the lengths and Ly determining the unit (smallest 
periodic) cell satisfy Ly = ^/3Lx while Ly = Lx/V^ for 
the class 2 . Accordingly, the Brillouin zone and Dirac 
cones (at the K-points) of these two lattice classes have 
different properties: Ky = Ky = 0 for the class 1 while 
Ky = —Ky = 27rl3Ly for the class 2 as schematized as in 
Fig. 3. Similarly to the case of commensurate lattices, if 
the channel is made of two layers of the same class, the 
Dirac cones of two layers are located at the same k-point 
and the transport gap is zero (not shown). However, if 
the channel is made of layers of different classes, there is 
a misalignment of Dirac cones of two layers in the k-space 
and the transport gap is finite. This essentially explains 
our obtained results. 

In the table I, we show the detailed description of some 
concrete devices investigated in this work and their trans¬ 
port gap. Ly is the average periodic length of two lay¬ 
ers along the Oy axis and 77 denotes the mismatch be¬ 
tween Ly and i.e., r] = \Ly — Ly\/Ly. As mentioned 
above, our calculations were performed with r] < 2% in 
all cases. First, we find that there are various possibili¬ 
ties of finite band-gap using this design strategy and Eg 



FIG. 4: Transport gap as a function of periodic length Ly 
along the Oy axis. 
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FIG. 5: Strain effects on the transport gap in different devices. 
The strain of angle 0 with respect to the transport direction 
(Ox) is applied [17]. 

can reach a few hundreds meV. Important, the depen¬ 
dence of Eg on the twist angle is quite complicated but 
Eg is inversely proportional to Ly. This is essentially 
explained as follows. The value of Eg is basically pro¬ 
portional to the distance KKy between Dirac cones of 
two layers along the /c^-axis [la ns nj]. In the consid¬ 
ered channels, KKy is actually equal to 27r/3I/^, which 
fully explain the property above. On this basis, we can 
also analytically compute the value of Eg. Note that 
in the low energy regime (i.e., \E\ 1 eV around the 

Dirac point), the energy dispersion of graphene is well 
described by the simple formula E{^ = Ehvpk with 
vf = 3aoto/2h and k = k — K (i.e., to = 2.7 eV and uq 
= 0.142 nm). The bandegde is thus defined as a func¬ 
tion of ky by Eedge{ky) = E3toao\ky\/2. In the consid¬ 
ered channel, the Dirac cones of two layers are located at 
different /c^-point, i.e., at Ky = 0 and Ky = 27r/3L^ re¬ 
spectively, and the energy-gap is thus determined as the 
value of Eedge{ky) at the middle point [HI |T 6 l [II]: he., 
ky = iT/3Ly (see also in Figs. 2(b,c,d)): 

Eg = ( 1 ) 

Ly 

In Fig. 4, we summarize the data obtained in all studied 
devices and, indeed, it confirms that the results match 
very well with the analytical formula above. 

Now, new questions arise. Since the dependence of Eg 
on the twist angle is quite complicated, it shall be diffi¬ 
cult to control this energy-gap in practice. Additionally, 
it is also difficult to control well the commensurateness of 
two layers. To partly solve these issues, we propose to use 
strain engineering as an additional ingredient. As sug¬ 
gested in HZ], the strain effects are very efficient to open 
a finite energy-gap in devices made of commensurate lat¬ 
tices. Similarly, this technique can be used to open an 
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energy-gap in the channels where two incommensurate 
layers are in the same lattice class, i.e., when Eg = ^ 
without strain. Additionally, in the cases of device with 
a small gap (i.e., if the two layers are from different lattice 
classes and Ly is large), strain can be used to further en¬ 
large Eg. These points are indeed shown clearly in Fig. 5 
for three typical cases: commensurate (AM 7 ^ 4 @AM 4 j), 
incommensurate (AMi 4 ^i 4 @AM 9 ^i 6 ) with zero gap and 
incommensurate (AMii^ 20 @AMi 3 ^ 2 o) with a small gap 
without strain. The gap is significantly broadened, i.e., 
by a few hundreds meV, when a small strain of only 4 % 
is applied. 

In conclusion, we have investigated the transport prop¬ 
erties of vertical devices made of two incommensurately 
misoriented graphene layers. For a given Oxy coordinates 
with transport direction parallel to the Ox axis, there 
are two different classes of graphene lattice depending on 
their lattice symmetry and on the properties of their Bril- 
louin zone. In particular, the two Dirac cones are either 


located at the same ky-point {Ky = Ky = 0) or at dif¬ 
ferent ky-points {Ky = —Ky = 27 rjSLy). On this basis, 
a misalignment of Dirac cones of the two layers occurs 
and, hence, a significant energy-gap (i.e., a few hundreds 
of meV) of transmission is achieved in devices made of 
two layers of different lattice classes. The gap is shown to 
be inversely proportional to the periodic length Ly. We 
have also shown the possibility of using strain to enlarge 
the energy-gap opening in this type of device. Thus, our 
study suggest an alternative strategy to open an energy 
gap in graphene channels (without altering the graphene 
lattices) by controlling the misorientation of graphene 
layers, which should be helpful for broadening the prac¬ 
tical applications of graphene. 

Acknowledgment. This research in Hanoi is funded 
by Vietnam’s National Foundation for Science and Tech¬ 
nology Development (NAFOSTED) under grant number 
103.01-2014.24. 


[1] A. C. Ferrari et a/.. Nanoscale 7, 4598 (2015). 

[2] I. Meric, M. Y. Han, A. F. Young, B. Ozyilmaz, P. Kim, 
and K. L. Shepard, Nat. Nanotechnol. 3, 654659 (2008). 

[3] A. Alarcon, V. Hung Nguyen, S. Berrada, D. Querlioz, J. 
Saint-Martin, A. Bournel, and P. Dollfus, IEEE Trans. 
Electron Dev. 60, 985-991 (2013) 

[4] M. Y. Han, B. Ozyilmaz, Y. Zhang, and P. Kim, Phys. 
Rev. Lett. 98, 206805 (2007). 

[5] N. Kharche and S. K. Nayak, Nano Lett. 11, 5274 (2011); 
S. Tang et al, Sci. Rep. 3, 2666 (2013). 

[6] A. Lherbier et al, Nano Lett. 13, 1446 (2013); A. Zabet- 
Khosousi et al., J. Am. Chem. Soc. 136, 1391 (2014). 

[7] Y. Zhang, T.-T. Tang, C. Girit, Z. Hao, M. C. Martin, A. 
Zettl, M. F. Crommie, Y. R. Shen, and F. Wang, Nature 
459, 820 (2009). 

[8] J. Bai, X. Zhong, S. Jiang, Y. Huang, and X. Duan, Nat. 
Nanotechnol. 5, 190 (2010). 

[9] G. Fiori, A. Betti, S. Bruzzone, and G. lannaccone, AGS 
Nano 6, 2642 (2012). 

[10] L. Britnell et al, Science 335, 947 (2012); L. Britnell et 
al., Nat. Gommun. 4, 1794 (2013); W. J. Yu et al., Nat. 
Nanotechnol. 8, 952 (2013). 

[11] O. V. Yazyev and S. G. Louie, Nat. Mater. 9, 806809 

( 2010 ). 

[12] S. B. Kumar and J. Guo, Nano Lett. Nano Lett. 12, 


13621366 (2012). 

[13] J. Zhang, J. Gao, L. Liu, and J. Zhao, J. Appl. Phys. 
112, 053713 (2012) 

[14] G. Fiori, S. Bruzzone, and G. lannaccone, IEEE Trans. 
Electron Dev. 60, 268-273 (2013). 

[15] V. Hung Nguyen, H. Viet Nguyen, J. Saint-Martin, and 
P Dollfus, Nanotechnol. 26, 165201 (2014). 

[16] M. Ghung Nguyen, V. Hung Nguyen, H. Viet Nguyen, 
and P Dollfus, Semicond. Sci. Technol. 29 (2014) 115024 

[17] V. Hung Nguyen, H. Viet Nguyen, J. Saint-Martin, and 
P Dollfus, Nanotechnol. 26, 115201 (2015). 

[18] A. Mishchenko et al., Nat. Nanotechnol. 9, 808813 
(2014). 

[19] G. R. Woods et al., Nat. Phys. 10, 451456 (2014). 

[20] V. M. Pereira, A. H. Gastro Neto, and N. M. R. Peres, 
Phys. Rev. B 80, 045401 (2009). 

[21] P. Moon and M. Koshino, Phys. Rev. B 85, 195458 

( 2012 ). 

[22] Z. F. Wang, F. Liu, and M. Y. Ghou, Nano Lett. 12, 
3833 (2012). 

[23] M. Koshino, New J. Phys. 17, 015014 (2015). 

[24] V. Hung Nguyen and P Dollfus, 2D Mater. (2015) in 
press; available at arXiv:1412.7583 



